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(57) Abstract 

Concentration of a target chemical, glucose, in the presence of another substance, NaCl, in a specimen (4) is determined by subjecting 
(2) the specimen (4) to radio frequencies (6, 16) up to about 5 GHz. The real and imaginary components of the reflected and/or transmitted 
signal are examined (18) to identify the presence and/or concentration of the chemical of interest The examination includes analysis of the 
effective complex impedance presented by the specimen (4) and/or the effective phase shift between the transmitted and reflected signals. 
The effects of NaCl on glucose concentration measurements can be nulled-out by examining impedance magnitude at a cross-over frequency 
or measuring NaCl concentration in a first frequency range and subtracting from a combined glucose/NaCl concentration measurement in 
a second frequency range. This technique can be used by diabetics to measure blood flucose in-vivo or in-vitro. 
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IMPROVING RADIO FREQUENCY SPECTRAL ANALYSIS 
FOR IN-VITRO OR IN-VIVO ENVIRONMENTS 

5 

FIELD OF THE INVENTION 
This invention relates generally to radio frequency spec- 
troscopy, and more particularly to improving specificity 
and accuracy of such analysis to determine the presence 
10 and/or concentration of a desired chemical among other 
substances within a specimen. 

BACKGROUND OF THE INVENTION 
Many conventional analysis techniques measure the concen- 

15 . tration of a chemical in a test specimen or sample, even 
where the specimen contains a complex mixture of chemi- 
cals. Such techniques include mass spectrophotometry, 
nuclear resonance, flame photometry, conductance and 
ref ractometry . While these techniques work, unfortunate- 

20 ly, their accuracy is tod often directly related to their 
cost. Further, many such techniques alter or destroy the 
specimen under test, and require relatively elaborate 
equipment . 

25 More recently attempts have been made to determine vari- 
ous properties of materials, using sound, electromagnetic 
waves, or single pulses as the basis for analysis. In 
contrast to conventional chemical analysis, wave and 
pulse-based techniques can provide a non- invasive in-vivo 

30 analysis. 
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For example, U.S. Patent no. 4,679,426 (July 1987) dis- 
closes a non- invasive in- vivo technique for measuring 
concentration of chemicals, sodium chloride for example. 
Periodic electromagnetic waves having a repetition rate 
5 of about 10 MHz to 100 MHz were coupled to a subject's 
finger, and sodium or chloride ions within the finger 
apparently distorted these waves. This distortion in the 
composite waveform was received from the finger, using 
the same electrode- antenna pair used to couple the waves 
10 to the finger. The composite waveform distortion was 

then examined, and found to provide meaningful data as to 
chemical concentrations. 

Glucose is an especially important chemical, a knowledge 
15 of whose absolute concentration level can be vital to 

diabetics. Several techniques for providing blood-sugar 
analysis are known, which permit subjects to determine 
their own glucose levels. Unfortunately many such tech- 
niques require invasive sampling of the subject. 

20 

One non- invasive technique for determining glucose levels 
in-vivo was disclosed in U.S. Patent no. 4,765,179 (Au- 
gust 1988) in which a periodic train of electromagnetic 
energy, preferably having a repetition rate of about 1 

25 MHz to 1 GHz, was coupled to a subject's finger. The 

composite waveform distortion was then analyzed and found 
to provide meaningful analysis of glucose levels in the 
range of about 50 to 150 mg percent. However, beyond 
about 110 mg percent, it was desirable to fine-tune the 

3 0 electromagnetic energy to maintain measurement accuracy. 
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Understandably, blood is a complex solution. Monitoring 
the concentration of glucose in blood presents substan- 
tial challenges to discriminate against other substances 
in the blood that may mask or alter the analysis results. 

5 

U.S. patent no. 5,508,203 described a non-invasive in- . 
vivo apparatus and method for determining a chemical 
level in a subject, including the chemical glucose. The 
use of frequencies up to about 1 GHz was disclosed and 
10 the disclosed apparatus permitted even lay persons in- . 

eluding diabetics to determine, for example, the level of 
glucose in their blood system. 

As useful as the invention disclosed in U.S. patent no. 

15 5,508,203 is, applicants have since realized that elec- 
trolytes, e.g., NaCl, KCi, Na 2 HP0 4 , and KH 2 P0 4 of varying 
concentrations in human blood can affect the accuracy of 
glucose measurements using that invention. In the human 
system, glucose concentrations typically range 60 mg/dl 

20 to about 150 mg/dl for a non-diabetic, and range from 
about 50 mg/dl to 500 mg/dl. In the human population, 
NaCl concentrations can range from about 135 mM to about 
145 mM. To effectively and confidently . measure glucose 
and/or its concentration in blood, a resolution of about 

25 10 mg/dl of glucose is desired. Non- invasive sophisti- 
cated laboratory grade test equipment can resolve glucose 
in-vitro to perhaps 1.5 mg/dl. Invasive consumer-grade 
can resolve glucose to perhaps 5 mg/dl with an accuracy 
of perhaps ±10%. Applicants are not aware of existing 

30 non- invasive devices for resolving glucose to the desired 
10 mg/dl level. 
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There is a need for a method and apparatus to reduce the 
varying concentration effects of electrolytes, especially 
NaCl, when measuring glucose concentrations in human 
blood. Such method and apparatus should be useable in- 
5 vitro and in- vivo, and should work in non- invasive in- 
vivo measurement environments. Further, such method and 
apparatus should be capable of use by lay persons. Such 
method and apparatus should also have applicability in 
measurements unrelated to analysis of bodily fluid, in- 
10 eluding applications in industry. 

The present invention discloses such a method and appara- 
tus . 

15 SUMMARY OF THE INVENTION 

A specimen containing a chemical of interest as well as 
other substances is via probes subjected to radio fre- 
quency electromagnetic signals having high frequency 
components extending to perhaps 5 GHz. Preferably such 

20 frequencies are sequentially presented using one sinewave 
frequency at a time, although simultaneously presented 
multiple frequencies may also be useful. Reflected 
and/or transmitted signal real and imaginary components 
at the specimen are then spectrally examined as a func- 

25 tion of frequency to identify the presence and/or con- 
centration of the chemical of interest. Such examination 
includes analysis of the effective complex impedance pre- 
sented by the specimen, and/or effective phase shift 
between the transmitted and reflected signal at the spec- 

3 0 imen. In this manner, greater specificity can be at- 
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tained with respect to detecting presence and/or concen- 
tration of a desired analyte or chemical of interest. 

For in-vitro measurements, a probe is inserted into the 
5 specimen and is coupled to a network analyzer, or similar, 
electronic system. In such in-vitro measurements, the • 
specimen may include blood or other bodily fluid, or may 
be a substance unrelated to bodily fluid. In in-vivo 
measurements, a network analyzer of similar electronic 
10 system may be coupled to electrode (s) on a probe. The 
probe is pressed against a subject ! s body, preferably a 
finger, and non-invasive analyses are made. 

Applicants have discovered that variable concentrations 

15 of electrolytes, especially NaCl, affect accuracy and 
specificity of glucose concentration measurements. At 
frequencies below about 1 GHz, increasing NaCl (or other 
small ions) concentration decreases impedance, whereas at 
higher frequencies the impedance is increased. Appli- 

2 0 cants believe that at the lower frequencies, ions can 

respond to the changing electromagnetic field adjacent to 
the probe ends, whereas this is more difficult at higher 
frequencies, whereat water dipoles appear to largely de- 
termine impedance. In general, applicants have learned 

25 that over a wide frequency regime, higher glucose concen- 
trations increases impedance, probably because the large 
glucose molecules hamper movement of electrolyte ions and 
water dipoles in a solution specimen. Of special inter- 
est, applicants have discovered that increasing NaCl 

30 concentrations over a wide frequency regime increase 
phase shift in a linear fashion, which phase shift is 
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insensitive to glucose concentrations. Using these dis- 
coveries, applicants can null -out or at least reduce or 
compensate for electrolyte concentration effects upon 
glucose concentration by using cross-over frequencies, 
and by examining different measurement parameters at 
different frequency regimes. 

In a blood specimen, electrolyte concentration effects 
are effectively "tuned out" by examining the magnitude of 
complex impedance using a cross-over frequency of approx- 
imately 2.5 GHz. This use of a cross-over frequency and 
complex impedance measurement provides low sensitivity to 
NaCl concentration and thus more accurate and specific 
glucose concentration readings. Such analysis improve- 
ment can be highly important, for example when the speci- 
men comes from a diabetic or suspected diabetic. 

Differential analyses may be made by combining impedance 
magnitude and phase shift measurement data. For example, 
high frequency phase shift measurements taken between 2 
GHz and perhaps 5 GHZ can provide data proportional to 
magnitude of ion concentration, particularly NaCl', On 
the other hand, impedance magnitude measurements made 
using lower frequencies, perhaps the 1 MHz to 400 MHz 
range, will provide a measure of combined concentration 
of glucose and ion concentration, again primarily NaCl. 
The high frequency phase shift data may be used to sub- 
tract out the effective NaCl concentration from the lower 
frequency impedance total concentration data. The result 
is a lower frequency measure of glucose concentration in 
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the specimen, a frequency regime in which measurement 
equipment is quite sensitive. 



Analysis equipment coupled to the impedance measurement 
5 data and phase shift measurement data can include look-up 
tables or the like, correlating phase shift data to NaCl 
concentration levels. For industrial applications, the 
look-up tables can store data correlating impedance, 
phase shift and frequency measurements to known substanc- 
10 es and concentration levels. This information can then 
be used to enhance nulling-out of NaCl in an impedance 
measurement made at a cross -over frequency. 

Output indicators coupled to such analysis equipment can 
15 enable even a lay user to readily understand what chemi- 
cal has been detected and at what concentration, or sim- 
ply to confirm that a safe concentration has been detect- 
ed for the chemical of interest. 

20 . Other features and advantages of the invention will ap- 
pear from the following description in which the 
preferred embodiments have been set forth in detail, in 
conjunction with the accompanying drawings. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a block diagram of a radio frequency spec- 
troscopy system; 

FIGURE 2 is a block diagram of the transmitter/receiver- 
30 signal processing system 14, shown in Figure 1; 
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FIGURE 3 A is a schematic of the calibration cell 66, 
depicted in Figure 2; 



FIGURE 3B is a Smith chart impedance versus frequency 
5 representation of the equivalent circuit depicted in 
Figure 3A; 

FIGURES 4A, 4B, and 4C depict signal amplitudes provided 
, by the system of Figure 1 for different target chemicals 
10 in analyte test solutions; 

FIGURE 5A depicts an in- vitro application of a radio fre- 
quency spectroscopy system with enhanced analysis sensi- 
tivity, according to the present invention; 

15 

FIGURE 5B depicts an in-vivo application of a radio fre- 
quency spectroscopy system with enhanced analysis sensi- 
tivity, according to the present invention; 

20 FIGURE 6A compares non-invasive and invasive impedance 

magnitude test data for a subject, using a test configu- 
ration according to Figure 5B; 

FIGURE 6B shows correction for electrolyte dilution for 
25 the sam data shown in Figure 6A; 

FIGURE 7A depicts linear relationship between electrolyte 
concentration and phase shift, independently of glucose 
concentration; 

30 
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FIGURE 7B depicts linear relationship between electrolyte 
concentration and phase shift in a PBS solution, 
independently of glucose and/or albumin concentration; 

5 FIGURE 7C demonstrates how improved specificity for a. 
target analyte can be realized by including measurements 
that are insensitive to a constituent in the specimen, 
for example, phase shift measurements at 1.5 GHz to null- 
out albumin concentration; 

10 

FIGURE 7D depicts a phase cross-over frequency of about 
20.1 MHz whereat phase shift data. is independent of glu- 
cose concentration; 

15 FIGURE 8A depicts the increase in impedance measured from 
about 0.1 MHz to about 1 GHz with increasing glucose 
concentration; 

FIGURE 8B depicts a frequency regime in which increasing 
2 0 NaCl and glucose concentrations increase impedance; 

FIGURE 8C depicts a frequency regime in which increasing 
NaCl concentration does not substantially affect 
impedance, but increasing glucose concentration increases 
25 impedance; 

FIGURE 8D depicts a 2.0 GHz to 2 . 1 GHz frequency regime . 
in which increasing NaCl concentration decreases imped- 
ance, while increasing glucose concentration increases 
30 impedance reasonably linearly; 
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FIGURE 8E depicts the non- linear behavior of impedance 
magnitude data over a 2.25 GHz to 2.75 GHz frequency 
regime as NaCl concentration is varied, according to the 
present invention; . 

5 

FIGURE 8F depicts the existence of a cross-over frequency 
at about 2.5 GHz at which NaCl concentration effects upon 
measured impedance are nulled-out; 

10 FIGURE 8G depicts frequency versus impedance changes for 
a specimen containing various substances, and demon- 
strates a possible gamma globulin saturation region. 



15 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Figure 1 depicts a radio frequency ("RF") spectroscopy 
system 2 for determining the presence of one or more 
target chemicals (depicted as x, y) in a cell membrane 
specimen 4, e.g., a human finger. The specimen finger 4 

20 is pressed against a probe pair 6, preferably disposed 

within a concave depression 8 formed in a lucite base 10. 
Probe pair 6 comprise two conductive rods that protrude 
slightly from the depression 8, permitting electrical 
contact to be made when finger 4 is pressed against the 

25 rods. Preferably the rods are brass, perhaps 0.2" (5 mm) 
outer diameter and protrude outward from the concave 
surface about 0.05" (1.3 mm), and into the lucite base 
about 0.5" (12 mm). Of course other tissue could be 
probed, e.g., an ear, and the specimen need not be a 

30 human. 
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A pair of transmission lines 12 electrically couples the 
electrode pair to a system 14 that includes a transmitter 
unit 16 and receiver-signal processor unit 18. Briefly, 
unit 16 transmits a high frequency signal via transmis- 
5 sion lines 12. to probes 6, which couple the signal to the 
specimen finger 4. Although the precise mechanism is not 
fully understood, it appears that the presence of target 
chemicals, e.g., x and/or y, within the specimen may 
cause energy transfer of certain spectra of the source 

10 signal from transmitter 16. The result is that a return 
signal from the specimen, present at pirobe pair 6 and 
coupled via transmission lines 12 to unit 18, differs 
from the source signal. Of course separate probe units 6 
could be used to couple the transmitter unit 16 to the 

15 specimen, and to couple the return signal from the speci- 
men to unit 18 . 

Unit 18 receives and processes the return signal such 
• that spectral signatures associated with the presence and 

20 concentration of various target chemicals within the 

specimen can be recognized. The processed data is then 
coupled to a display system 2 0 that conveys the detected 
information to a user. Operation of the receiver-signal 
processor unit. 18 can be tailored, manually or automati- 

25 cally by a neural network, to recognize specific target 
chemicals, for example glucose within the blood stream 
within- finger specimen 4. In such instance, the various 
output devices within display system 20 might provide a 
user with calibrated data as to his or her glucose con- 

30 centration. 
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Display system 20 may include a monitor that can display 
a spectrum analyzer output (22A) , and/or alpha- 
numeric/graphical output (22B) . Display system 2 0 may 
also include a bar graph or alpha-numeric indicator 24 
5 indicating, for example, the concentration level of the 
target chemical, for example, glucose. A calibrated 
output meter 2 6 could provide the user with concentration 
data. Alternatively, a simple "GO/NO GO" output indica- 
tor 28 could alert the user that excess glucose concen- 
10 tration has been detected. A diabetic user would thus be 
alerted to take insulin immediately. 

Figure 2 is a block diagram of the transmitter/receiver- 
signal processing system 14. Oscillator 50 generates a 

15 high frequency stimulus signal that will be transmitted 

via probes 6 to specimen 4. In the preferred embodiment, 
oscillator 50 provides a 30 MHz fundamental square wave 
having a 50% duty cycle, and transition times of a few 
nanoseconds. As such, the oscillator output frequency 

20 spectrum will be rich in harmonics, the odd-numbered 
harmonics predominating. In the frequency domain, a 
perfect square-wave source signal would have harmonics 
with a sin(x)/x envelope, where x represents a harmonic 
frequency. 

25 

The spectral output of such an oscillator 50 is commonly 
referred to as a comb spectra, as the various spectra are 
uniformly spaced similar to the teeth on a comb. The 
power output level at the oscillator output is preferably 
30 about 1 mW, which is 0 dBm, although other power levels 
may also be used. 
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In the preferred embodiment, the various source signal 
spectra are harmonically related since generation of a 
•pulse train provides the harmonic frequencies, automati- 
cally. However the source frequencies need not be har- 
5 monically related, and a single oscillator 50 may be 

rapidly changed between discrete frequencies (e.g., in • 
the manner of spread spectrum transmitters) . Alterna- 
tively, oscillator 50 could comprise a plurality of sig- 
nal generators whose separate frequency outputs may or 

10 may not be harmonically related. If harmonically relat- 
ed, one such generator could provide a sinusoidal output 
at a fundamental frequency, e.g., 30 MHz . A second gen- 
erator could provide a 60 MHz sinusoidal output, a third 
generator could provide a 90 MHz sinusoidal output, and 

15 so forth. In a different embodiment, one such generator 
might provide an output at frequency fl, a second genera- 
tor might provide an output at f2, not harmonically re- 
lated to fl, and so forth. 

20 As used herein, oscillator 50 is understood to be a 

source of electromagnetic signal that contains a plurali- 
ty of high frequency components, regardless of whether 
such components represent harmonics of a single source 
frequency, or represent many source frequencies, that 

25 need not be harmonically spaced-apart . 

Unit 52 preferably includes an amplifier stage and a 
power splitter, and comprises a MAR-3 amplifier and a 
Cougar amplifier stage and a power splitter in the pre- 
30 f erred embodiment. These commercially available compo- 
nents boost the oscillator signal provided to divider 54 
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to about 15 dBm, and provided to power splitters 62, 64 
to about 3 dBm. In turn, each power splitter 62, 64 
divides the thus amplified signal into two signals at 
nodes A and B, each having 0 dBm power output. Splitters 
62, 64 are preferably wideband, e.g., about 10 MHz to 
1, 000 MHz (or 1 GHz) . 

The intermediate frequency ("IF") for system 14 is pref- 
erably 21.4 MHz, an intermediate frequency commonly used 
in commercial equipment, for which frequency many stan- 
dardized transformers and circuits are readily available. 
High- side mixing injection preferably is used. Thus, to 
generate a local oscillator frequency that is 21.4 MHz 
higher than a center frequency, it is necessary to devel- 
15 . op a synthesized reference 6.4 MHz signal. Unit 54 di- 
vides the fundamental frequency of the oscillator signal 
by 6, to yield a nominal 5.0 MHz reference signal. 
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This 5.0 MHz reference signal and a 6.4 MHz phase-locked 
crystal controlled oscillator signal 58 are processed by 
offset module 56. : Offset module 56 outputs on line 60 a 
signal having a frequency of 6.4 MHz that is phase locked 
to the 30 MHz frequency of oscillator 50. Because phase 
lock loop systems are well known in the art of digital . 
25 signal processing design, further details of the: genera- 
tion of the frequency locked 6.4 MHz signal on line 60 
are not presented here. 



30 



In Figure 2, calibrator unit 66 is an electronic model of 
a typical human finger, essentially the electronic equiv- 
alent circuit of a finger specimen 4. While calibration 
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unit 66 approximates the specimen impedance, unit 66 will 
not include the target chemical. 

Figure 3A details the circuitry within calibration unit 
5 66, namely two segments of transmission line having 50 Q 
impedance at 400 MHz, and assorted resistors and capacir 
tors. The transmission lines, . resistors and capacitors 
were selected empirically by comparing frequency versus 
impedance data from human fingers with data from the 

10 equivalent circuit of Figure 3A. Figure 3B is a Smith 
chart impedance versus frequency representation of the 
equivalent circuit of Figure 3A. Point A in Figure 3B 
represents an impedance of about 192 Q/-201 Q at 10 MHz, 
point B is 39.5 Q/11.5 Q at 300 MHz, C is 52 Q at 400 

15 MHz, and point D is about 57 Q/-2.6 Q at 500 MHz. 

With further reference to Figure 2, as will now be de- 
scribed, various components are replicated to provide a 
processing path for the transmitted source signal, and to 

20 provide a processing path for what will be termed the 

sampled return (or received) signal. The sampled return 
signal advantageously permits compensating the system of 
. Figure 2 for component variations and drift between what * 
will be termed the received and the transmitted signal 

25 processing paths. 

More specifically, the response of specimen 4 to the 
source signal (e.g., the return signal at the probe pair 
6) is switchably sampled by switch SI with the response 
30 of the calibration unit 66 to the source signal. Harmon- 
ic frequency-by- frequency, the output from probe pair 6 
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and from calibration unit 66 are sampled, the output of 
SI providing a sampled return signal at node C to the 
remainder of system 18, Of course, if source oscillator 

50 provided discrete frequencies that were not harmon- 
5 ically related, it is understood that f requency-by- fre- 
quency, the output from probe pair 6 and from calibration 
unit 66 would be sampled. In, the preferred embodiment, 
the frequency bands of interest begin with about the 
sixth or seventh harmonic of source oscillator 50, e.g., 

10 about 195 MHz, and extend to about 1 GHz, or higher, 

which range is the bandwidth of system 18. Within that 
bandwidth, individual frequencies are sampled between 
probe pair 6 and calibration unit 66. 

15 Switch SI preferably is a commercially available mono- 
lithic microwave integrated circuit ("MMIC"), a relay, or 
other switching mechanism. SI switches between the probe 
6 output and the calibrator under control of a micropro- 
cessor 74 within system 14. In the preferred embodiment, 

20 microprocessor 74 was a Motorola 68HC11, although other 
microprocessors could be used instead. 

51 may sample the output of probe 6 for a time period 
ranging from perhaps 30 ms to perhaps 7 seconds, and then 

25 may sample the output of the calibration unit 66 for a 

time period also within that range, the duty cycle typi- 
cally being aperiodic. For example, during the time SI 
is coupled to probe 6, the probe output signal is sampled 
for one or more frequencies that are harmonics of the 

30 fundamental frequency of oscillator 50 (or for one or 
more discrete frequencies provided by an oscillator 50 
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that does not provide harmonics) . During the time SI is 
coupled to the calibration unit 66, the response of cali- 
bration unit 66 to one or more frequencies that are har- 
monics of the fundamental oscillator 50 frequency are 
5 sampled. 

Understandably, if components 76T and 76R, 78T and 78R, 
80T, 80R, 90T and 90R (to be described) were identical 
and exhibited no drift, calibration unit 66 could be. 

10 dispensed with, and SI replaced by a wire making a perma- 
nent connection in the probe 6 SI position. Such an 
ideal system would require no mechanism for compensating 
for drift and other differences in the signal processing 
paths for the harmonics of the oscillator signal 50, and 

15 for the harmonics in the return signal obtained from . 
probe 6 . 

In practice, variations in temperature and/or pressure 
between probe pair 6 and the tissue in the specimen 4 may 

20 contribute some error to the measurement process. To 

permit microprocessor 74 to compensate for such error, in 
addition to providing the microprocessor with phase and 
amplitude information for harmonics, phase and amplitude 
information is also provided for the oscillator funda- 

25 mental frequency. This frequency has been found experi- 
mentally to be sensitive to such temperature and/or pres- 
sure variations. It is understood that suitable tempera- 
ture and/or pressure transducers and analog- to-digital 
conversion components that are not shown in Figure 1 are 

30 used. 
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As shown in Figure 2, within the transmitted source sig- 
nal processing path, a bandpass filter 68T has a center 
frequency equal to that of oscillator 50, e.g., 30 MHz, 
and a bandwidth of about 1 KHz to perhaps 1 MHz . Other 
bandwidths could be used and in fact, a 30 MHz lowpass 
filter might instead be used. The transmitted signal . 
from node A is coupled to bandpass filter 68T, and the 30 
MHz center frequency component of this signal passes from 
filter 68T and is amplitude limited by limiter 70T. The 
thus bandpass filtered and amplitude limited signal is 
coupled to an input of a phase detector 72 . 

In a parallel path, the sampled return signal from switch 
SI, present at node C, passes through a similar 30 MHz 
bandpass filter 68R, amplitude limiter 70R to provide a 
second input to phase detector 72. (The letter T or R 
attached to a reference element herein denotes that the 
element is used in the transmitted source path, e.g., 
68T, or is used in the sampled return signal path, e.g., 
68R.) 

Phase detector 72 compares the difference in. phase be-'", 
tween the transmitted 30 MHz fundamental frequency. and 
the sampled return 30 MHz fundamental frequency signal. 
The phase detector 72 output signal voltage will be pro- 
portional to such phase shift, e.g., a number of mV per 
each degree of phase shift. As shown in Figure 2, the 
phase output information from detector 72 is coupled to 
microprocessor 74 for analysis. 
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Proceeding horizontally across the top of Figure 2, par- 
allel paths are also depicted for processing the trans-, 
mitted source signal harmonics (available at node B) and 
the sampled return signal harmonics from switch SI 
5 (available at node C) . These two horizontal paths use 
substantially identical components (as denoted by the 
nomenclature) to provide transmitted and sampled return 
signals at an intermediate frequency (IF) that is about 
21.4 MHz in the preferred embodiment. 

10. 

Briefly, the components now to be described resolve the 
harmonic frequency components of the signals at node B 
and node C into preferably four bands of discrete fre- 
quencies, depending upon what harmonics of the source 
15 oscillator signals are desired to be examined. 

Much of the remainder of the signal processor functions 
as a scanner-receiver, that under microprocessor control 
scans discrete harmonic frequencies of interest. The 

2 0 transmitted source signal path components will first be 
described, it being understood that identical components 
are used in the parallel sampled return signal path, as . 
indicated by the nomenclature, e.g., 76T, 76R, 78T, 78R, 
etc. Bandpass filter 76T (and thus also 76R) preferably 

25 ; is a filter bank that includes an internal MMIC switching 
mechanism operating under control of microprocessor 74. 
The input port of filter 76T passes the transmitted sig- 
nal from node B through an internal switch into two banks 
of pre-shaping three-pole bandpass filters. These first 

30 two internal filter banks have bandpasses of 195 MHz to 
395 MHz, and 395 MHz to 805 MHz. Still within filter 
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bank 76T, the outputs from the 195-395 MHz and 395-805 
MHz filters pass through additional internal MMIC switch- 
es and bandpass filters. These additional filters pass 
195-295 MHz , 295-405 MHz, 405-610 MHz, and 605-815 MHz . 
Still within 76T, the variously filtered components are 
combined into a single signal that is amplified by ampli- 
fier 78T. 

In similar fashion, the sampled return signal at node C 
is passed through switching bandpass filters within band- 
pass filter bank 76R, and the variously filtered compo- 
nents are combined and amplified by amplifier 78R. While 
the operation of bandpass filter banks 76T, 76R has been 
described with reference to specific frequency bands, 
those skilled in the art will recognize that the frequen- 
cies comprising the signals at nodes B and C may be fil- 
tered using bandpass filters having different ranges of 
bandpass. Because the design of units 76T, 76R is known 
to those skilled in the relevant art, schematics are not 
here provided. 

For example, if the target chemical of interest is best 
resolved by examining say the seventh harmonic of the 30 
MHz transmitted source signal (or a given discrete fre- 
quency of a source signal providing a plurality of fre- 
quencies not necessarily harmonically related) , micropro- 
cessor 74 is caused to control the switching within units 
76T, 76R to pass 210 MHz frequency components, e.g., to 
select the 195 MHz-295 MHz bandpass. Amplifiers 78T, 78R 
preferably have sufficient gain to compensate for attenu- 
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ation caused by filters 76T, 76R, and have a bandwidth of 
at least 195 MHz to 815 MHz. 



Of course, if amplifiers 78T, 78R were ideal and not 
5 subject to front -end overload, it would be possible to 
delete the bandpass filter systems 76T, 76R, and rely 
upon the operation of mixers 80T, 80R, and narrow band IF 
units 90T, 90R (to be described), to separate the various 
harmonic components of the oscillator signal and of the 
10. return signal. . 

As shown by Figure 2, the output signals from amplifiers 
78T, 78R are provided as an input signal to mixers 80T, 
8 OR. Frequency synthesized local oscillators LOl or L02 
15 provide respective second input signals to mixers 80T, 
80R, via a.MMIC switch S2 (or similar device) that 
. switches between the two synthesized oscillator signals 
under control of microprocessor 74 . 

20 The synthesized LOl or L02 signals are then frequency 
mixed against the selective spectral components of the 
transmitted source signal and sampled return signal that 
have been switchably selected to pass through filter 
banks 76T, 76R. The LOl or L02 output signals are 21.4 

25 MHz above the harmonic frequency of interest. Because of. 
the difficulty associated with implementing a synthesized 
local oscillator whose output frequency can range from, 
about 231.4 MHZ (e.g., 7x3 0 MHz + 21.4 MHz) to perhaps 
800 Mhz (e.g., about the twenty-sixth harmonic 26x30 MHz 

30 +21.4 MHz), the preferred embodiment employed two local 
oscillators, LOl, L02 . If, however a suitable synthe- 
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sized oscillator having a two-octave frequency output 
could be implemented, such oscillator would replace LOl, 
L02 and the necessity for S2 . 

Stages 90T, 90R are narrowband intermediate frequency 
circuits that pass a 21.4 MHz center frequency with a 
bandwidth of about 25 KHz . Of course by suitably offset- 
ting mixing frequencies, an IF of other than 21.4 MHz 
could be used. In the preferred embodiment, IP units 
90T, 90R are similar to IF units commonly found in com- 
mercially available cellular telephones. 

The harmonic frequency information passing through IF 
units 90T and 90R are input to phase detector 92 . Phase 
detector 92 compares transmitted source and sampled re- 
turn signals at each harmonic frequency of interest. 
The difference in phase between these signals is then 
provided by phase detector 92 to microprocessor 74 . At 
the same time, the relative voltage levels from the IF 
units 90T, 90R at node D are also provided (after suit- 
able analog to digital conversion, converter not shown) 
to microprocessor 74. 

To recapitulate, microprocessor 74 receives phase infor- 
mation from detector 92 that is relative to the various 
harmonics of the source signal (or discrete frequencies 
of interest if a non-harmonic generator 50 is employed), 
and that is relative to the various harmonics (or dis- 
crete frequencies) of the source signal as altered by the 
target substance and received at the probe pair 6. Simi- 
larly, microprocessor 74 receives amplitude information 
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of IF units 90T and 90R relative to the various harmonics 
(or" discrete frequencies of interest) of the source sig- 
nal, and that is relative to the various harmonics (or 
discrete frequencies) of the source signal as altered by 
5 the target substance and received at probe pair 6 . Fur- 
ther, to permit compensation for probe temperature and/or 
probe -specimen pressure variations, limiters 70T, 70R 
provide microprocessor 74 with amplitude of the source 
frequency, and with amplitude of the source frequency as 
10 altered by the target substance and received at probe 

pair 6, while detector 72 provides similar phase infor- 
mation for the source frequency. 

Microprocessor 74 operates under program control, gener- 
is ating data for further processing by a so-called neural 
network, look-up table, algorithm, or other method of 
signal processing, shown symbolically in Figure 2 as 
element 100. In a manner known to those skilled in the 
relevant art, a neural network 100 can be "trained" to 
20 recognize a spectral signature associated with a given 

target chemical, glucose for example. To ease this rec- 
ognition, neural network 100 can optimize the manner of 
signal processing within unit 14. 

25 For example, the operation of filter banks 76T, 76R can 

be altered under control of microprocessor 74 . In a more 
generalized embodiment, the number and bandwidth of indi- 
vidual bandpass filters within units 76T, 76R could be 
dynamically modified by suitable MMIC-selection, all 

30 under microprocessor control. However, unit 100 may 

simply be a look-up table, correlating relative amplitude 



23 



WO 00/09996 PCT/US99/18096 
changes in a return signal with harmonic frequency 
against presence or concentration of a target chemical in 
the specimen. Further, a suitable neural network 100 
might control microprocessor 74 to optimize the genera- 
tion of discrete frequencies, based upon processed signa- 
ture data. For example, if a certain set of frequencies 
from oscillator 50 provided a slight spectral signature, 
network 100 might direct oscillator 50 to provide slight- 
ly different frequencies until the signature was more 
recognizable. 

Microprocessor 74 in turn provides output signals to 
output indicator (s) 20. As has been described, output 
indicator (s) 20 can, in a variety of formats, display 
information enabling a user to determine the presence and 
concentration of a desired target chemical (e.g., x) in a 
specimen. In the preferred embodiment, the specimen is 
in fact a finger of the individual using the disclosed 
system. Although the system shown in Figures 1 and 2 was 
implemented in breadboard fashion, those skilled in the 
art will appreciate that it may in fact be fabricated in 
a handheld, battery operated, portable unit. In such em- 
bodiment, output indicator (s) 20 would preferably include 
liquid crystal displays (LCDs) or simple GO /NO GO indica- 
tors, to preserve power and space. Preferably base 10 
would be attached to the case housing the remainder of 
the system for ease of portability. 

Figures 4A and 4B represent multiple averaged in-vitro 
data obtained with the system of Figures 1 and 2, using 
as a test specimen whole blood (e.g., red blood cells) to 
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which glucose or lactose or sucrose or urea or NaCl was 
added as a test chemical . The test cells were compared 
to a calibrated cell that contained only red blood cells. 
Figure 4C represents similar data for whole sheep's blood 
5 (e.g., no glucose) , and for sheep's blood with various . 
concentrations of glucose, where the nomenclature "Blood 
102" denotes 102 mg-% or 102 mg per dL glucose. Typical- 
ly, a healthy human has perhaps 80-120 mg% glucose, while 
a diabetic has 200-400 mg-% glucose. The vertical axis 
10 in Figure 4C represents the vector amplitude the return 
signal, taking into account magnitude and phase. The 
horizontal axis represents harmonics of a 30 MHz source 
frequency, the first harmonic being at 210 MHz. 

15 . To minimize probe-related variables, the specimens in 
Figures 4A, 4B and 4C were tested using parallel . plate 
capacitive cells. These cells comprised two dielectric 
substrates having a relative permittivity approximating 
that of water (= 80) , with an electrode surface baked 

2 0 onto each substrate. The test substance was placed in a 
chamber between the substrates. 

The varying degree of signal amplitude shown in Figures 
4A, 4B, and 4B are termed "spectral signatures". What is 

25 depicted is the difference in amplitude between the cali- 
. brated cell (analogous to the use of the calibration unit 
66 in Figure 2) and the test specimen (analogous to the 
use of probes 6 and specimen 4 in Figure 1) . These data 
indicate that the system of Figures 1 and 2 may be used 

30 to discern the presence of a target chemical within a 
test specimen or sample. 
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A preferred application is the detection of excess glu- 
cose in a user's blood, e.g., within the specimen. Be- 
cause the present invention operates non-invasively, it 
suffices for the user to press his or her finger against 
the probe pair 6, as shown in Figure 1. In response to 
the high frequency, high harmonic content signal from 
transmitter 16, chemicals within the specimen can recog- 
nizably cause energy transfer of certain spectral compo- 
nents of the transmitted source signal. It is hypothe- 
sized that within the specimen, the target chemical glu- 
cose interacts with the lipid bilayer and/or red blood 
cell membranes. 

Thus, in the presence of frequency components from the 
signal transmitted via probes 6, the glucose seems to 
bring about non-linear intermodulation or mixing of fre- 
quency components, possibly due to a non- linear dielec- 
tric phenomenon involving capacitance associated with 
glucose. Using the system of Figures 1 and 2, a diabetic 
user may rapidly obtain glucose concentration level in- 
formation. Signal processing by unit 18 would, essen- 
tially in real time, provide glucose level information on 
display unit 20 . 

Of course other target chemicals may also be detected, 
including for example fructose, galactose, alcohol. For 
example, a system according to what is disclosed herein 
may be used to sense alcohol in a motorist's system, 
either by a motorist before attempting to drive, or by a 
police officer attempting to determine whether an indi- 
vidual is under the influence of alcohol. 
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Because the disclosed system of Figures 1 and 2 appears 
to. be sensitive to boundary conditions at a lipid bilayer 
membrane, disruptions to such boundary conditions may be 
detected by a spectral signature. Thus, the presence of 
glucose in varying amounts at a membrane may be detected. 
In a different utility, however, trauma to a specimen 
that interferes with such boundary conditions may also be 
detected, primarily for the purpose of providing medical 
treatment. For example victims of electrocution may 
received localized injuries, for example on an arm. 
Unless the injury sites are promptly treated by the in- 
jection of certain medication that is potentially rather 
toxic, the victim will lose the injured limb or die. Use 
of the invention disclosed herein would permit diagnosis 
of such injury sites, and quantizing the injury to facil- 
itate prompt and accurate medical treatment. 

Subsequent to the invention described with reference to 
Figures 1-4B, applicants came to appreciate the role that 
changing electrolyte concentrations can have upon glucose 
concentration measurements in blood specimens. Appli- 
cants further discovered that it is possible to improve 
analysis for a desired chemical by reducing the effects 
upon such analysis of varying concentrations of other 
substances in the specimen. 

Figures 5A and 5B respectively show in-vitro and in-vivo 
applications of improved analysis using a system 200, ac- 
cording to the present invention. In Figure 5A, pref- 
erably two probes 202A, 202B are coupled by short lengths 
of coaxial cable 12 to ports A and B of a frequency gen- 
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erator and analyzer system 250. In general, the trans- 
mitted signal is sent from port A or port B, and a por- 
tion of the transmitted signal is reflected by the speci- 
men back into the transmitting port. In transmission 
5 mode (e.g., Figure 5B) , port B returns the fraction of 
the signal transmitted via port A through the subject's, 
finger. 

In the embodiments of Figures 5A and 5B, cables 12 pref- 
10 erably are 20 cm or less lengths of coaxial cable, and 
probes 2 02A, 202B preferably are Hewlett Packard HP 
85075B dielectric probes. These probes are coaxial in 
construction, having an outer diameter of perhaps 2 cm 
and a probe length of perhaps 3.8 cm. The probes have a 
15 center conductor that is surrounded by a groundplane 

sheath at the probe tip. However, other cable couplings 
and probes could also be used. 

As will be described, system 250 includes a transmitter 
20 unit 260 that can output discrete sinusoidal waveforms 
that are spaced-apart in frequencies linearly or 
logarithmically in user-selectable steps. Further, the 
output frequencies are stepped between user-selectable 
lowermost and uppermost frequencies f 1 and f u , respec- 
25 tively. In the preferred embodiments, t x was about 300 

KHz, f u was about 3 GHz, with approximately 801 linearly- 
spaced frequencies output between f x and f u . Applicants 
believe, however, that an f u of about 5 GHz would also be 
useful to the present invention. In the preferred em- 
30 bodiment, system 250 was implemented using a commercially 
available Hewlett Packard HP 8753A network analyzer with 



28 



WO 00/09996 PCT/US99/18096 
an HP 85046A S-parameter test set. However, other sys- 
tems implementing similar functions could* be used in- 
stead. 

5 System 270 further contains a receiver and signal proces- 
sor unit 270 that analyzes waveforms associated with 
signals transmitted by and/or at least partially reflect- 
ed back to system 270. The waveforms under analysis are 
associated with discrete user-programmable frequencies. 

10 The analysis can examine real and imaginary components of 
these waveforms, including complex (e.g., having real and 
imaginary components) reflection coefficient data. These 
various data are. signal processed by unit 270 to provide 
information including complex impedance magnitude (Z) , 

15 phase shift, and/or permittivity. 

Among the electrolytes, NaCl has the most significant 
influence on measurements, in that its normal concentra- 
tion range in the human body is 135-145 mM (millimolar) , 

20 whereas KC1 , by example, is only abut 4-10 mM. Substanc- 
es such as urea were confirmed to not influence glucose 
measurements, probably because urea has a molecular size 
that is one-third that of glucose, and has a 
physiologically controlled concentration ranging from 5- 

25 40 mg/dl. The range of glucose in a human normally is 

about 50 mg/dl (or mg%) to 150 mg/dl, and can reach about 
500 mg/dl in a diabetic. 

In Figure 5A, probe 202 contacts a specimen of interest 
30 204, perhaps about 40 ml, retained within a beaker or 

receptacle 206 whose volume is perhaps 100 ml. Specimen 
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204 includes a chemical of interest denoted X, as well as 
one or more other substances, denoted collectively Y. In 
a preferred embodiment, specimen 204 is a bodily fluid, 
for example blood, X is glucose (whose presence and/or 
5 concentration is to be determined) , and Y may include 
varying concentrations of blood electrolytes such as 
NaCl, Na 2 HP0 4 , KCl , and KH 2 P0 4 , as well as proteins and 
lipids. 

10 Although large concentrations of proteins and lipids are 
also found in blood, the human body maintains relatively 
tight control over variations in such substances, and 
thus their presence appears not to substantially affect 
measurements according to the present invention. 

15 

In an industrial application, specimen 204 may be a solu- 
tion in which X and Y represent different chemicals, in 
which the presence and/or concentration of X is to be 
20 discerned, for example to confirm quality control of the 
production of solution 204. 

A second container 210 into which probe 202B is inserted 
contains a test or control solution 208 that intentional - 
25 ly lacks at least one chemical found in specimen 204. 

Both specimens preferably are retained at a same tempera- 
ture by partially immersing containers 206, 210 in a 
preferably constant temperature bath 212 maintained with- 
in a larger beaker or container 214. 

30 
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In Figure 5A, analyzer unit 250 is operated with signals 
at ports A and B in a reflectance mode, e.g., in which 
signals transmitted out of each port are at least par- 
tially reflected back into the ports by the respective 
5 . specimens. From the real and imaginary components of the 
reflected signal data, useful information as to the pres- 
ence and concentration of at least one chemical in solu- 
tion 204 may be . determined, according to the present 
invention. 

10 ' 

Applicants have discovered that the real and imaginary 
components of the reflected signals can be affected by 
the nature and content of the specimen solutions in the 
immediate vicinity of the tips of the probes. What is 

15 believed to occur is that fringing fields extend from the. 
center conductor of the preferably dielectric probes to 
the surrounding ground plane. As the properties of the 
specimen solutions change, e.g., due to the presence and 
concentration of one or more chemicals or other substanc- 

20 es therein, the fringing field is affected. The alter- 
ations to the fringing field in turn affect .the reflected 
signals being returned to ports A and/or B of the analyz- 
er unit 250. 

25 The complex data gathered and processed by unit 250 is 

coupled as input to a computer unit 2 80 for further pro- 
cessing. If desired, computer unit 280 may include any 
or all of the output indicators 22A, 22B, 24, 26, 28 
described earlier with respect to Figure 1, as well as 

30 any other output indicator (s) that may be desired. 
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Computer unit 2 80 may be a personal computer executing a 
software routine permitting conversion of the real and 
imaginary data it receives into forms including the mag- 
nitude of the effective complex impedance Z presented by 
the specimen, phase shift between signals transmitted and 
at least partially reflected back by the specimen, ef- 
fective permittivity, and the like. 

In the preferred embodiment, computer 280 executed Excel 
spreadsheet software to convert the incoming complex data 
into more useful form. A modified Bao procedure was 
adopted, in which complex impedance (Z) is determined 
from the complex reflection coefficient (r) at the inter- 
face between the flat end of a probe, e.g., 202A, and a 
specimen solution, e.g., 204. 

1 + r 

The characteristic impedance Z Q of coaxial line 12 may be 
calculated from the relationship: 



Z =377 
° \ 



e R 2*n 



(2) 



in which 377 represents impedance of air, Jb is the out- 
side diameter of the probe, a is the diameter of the 
inner lead on the probe, fi R is the permeability of air, 
and e R is the permittivity of Teflon. 



However, measured reflection coefficient from analyzer 
250 is not necessarily an accurate representation of r, 
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due to errors caused by the container 2 06, the coaxial 
line 12, and connectors at port A, for example. The Bao 
procedure reduces these errors, using a calibration pro-- 
cedure based on a linear assumption. This assumption and 
5 the values collected from the calibration procedure give 
rise to a matrix derivation 

e=-^-= — - (3) 

in which A* is a frequency dependent complex constant 
related to a scattering matrix. 

During the course of experimentation, applicants realized 
that if analyzer 250 were calibrated with port connectors 
and coaxial cables 12 attached, the analyzer output would 
be r, whereupon use of the Bao matrix procedure would be 
15 unnecessary. 

Thus, while equation (1) is valid, its real and imaginary 
components should be separated to be effectively used by 
computer 280 during execution of a data processing rou- 
20 tine, e.g., an Excel spreadsheet program. 

Consider then equations (4) and (5) , in which p is the 
complex reflection coefficient output from analyzer 250: 

z *e*i — 5 (4) 
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Imag , ? I 3 > 



Z Mag^ Z Lal^ 2 Imag (6) 



Euler's formula is used as shown in equations (6) and (7) 
to convert equations (4) and (5) to the more commonly^ 
encountered impedance magnitude and phase quantities: 

2 e = tan _1 — £££ (7) 

Z Raal 

Referring back to Figure 5A, the various analytes in a 
blood specimen, especially small ion electrolytes (also 
called blood electrolytes) , can measurably affect the im- 
pedance and phase angle. In an application in which 
glucose concentration is to be determined, what actually 
may be measured with system 2 00 is the effect of glucose, 
e.g., X in Figure 5A, upon ions or water dipoles in the 
specimen solution 204. Applicants have discovered at 
certain cross-over frequencies output by system 250, the 
effects of other substances Y in the specimen 204, in- - 
eluding electrolytes, can be reduced or nulled-out. For 
example, at a cross-over frequency of about 2.5 GHz, the 
concentration effects of NaCl and most probably other 
blood electrolytes in a blood specimen are nulled-out, 
without degrading glucose concentration measurements. In 
an analytical scheme in which N equations would have to 
be solved for N unknowns, the ability to null-out elec- 
trolyte concentrations effectively reduces the number of 
variables and thus the number of equations that must be 
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solved. The end result is that glucose concentration can 
be determined with higher specificity and confidence. 
Further, as described later herein, phase shift measure- 
ments (e.g., comparison between transmitted and reflected 
signals) over a wide frequency regime provide a surpris- 
ingly linear response to electrolyte concentration. The . 
phase shift data can then be used to compensate for NaCl. 
concentration contributions to total impedance measure- 
ments made at frequencies lower than the 2.5 GHz cross- 
over frequency. 

When generator 260 outputs frequencies greater than per- 
haps 1 GHz or so, the specimen impedance magnitude ap- 
pears to be primarily a function of the ability of water, 
dipoles to respond in the presence of the resultant os- 
cillating field in the vicinity of the probe (s). At 
output frequencies less than perhaps 500 MHz, impedance 
magnitude seems to be more a function of ionic response 
to the oscillating field in the vicinity of the probe (s) . 
Within a blood specimen, NaCl is an important source of 
such ions. At inbetween frequencies, the impedance func- 
tion transitions. 

Below approximately 500 MHz, glucose in the specimen 
solution appears to impede ionic mobility in responding 
to the oscillating field, and thus the effective imped- 
ance increases. For example, between about 10 MHz and 
100 MHz, impedance, change due to NaCl concentration 
changes in the specimen are substantially stronger than 
impedance changes due to concentration changes in glu- 
cose . 
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Applicants have discovered that. at test frequencies below 
about 1 GHz, increasing concentrations of NaCl decrease 
impedance magnitude ("Z"), and that at a cross-over fre- 
quency of about 2.5 GHz, impedance measurements are sen- 
sitive to glucose concentration but insensitive to elec- 
trolyte concentration. Further, applicants have learned 
that over a wide frequency regime, phase shift increases 
linearly with increasing NaCl concentration, with little 
or no effect due to changing glucose and/or albumin con- 
centration. Thus, it appears that at higher frequencies 
(e.g., above 1.5 GHz or so), larger molecules simply do 
not respond sufficiently rapidly to meaningfully influ- 
ence phase shift measurements. By contrast, 
electrolytes, including NaCl, have small ions that can 
respond measurably with respect to phase shift measure- 
ments. As described herein, collectively, these discov- 
eries provide measurement protocols to reliably and with 
specificity determine glucose concentration, despite the 
presence of electrolytes of varying concentrations. 

In Figure 5B, a non-invasive system for in-vivo testing . 
is depicted. In this embodiment, network analyzer or 
system 250, and computer system 280 may be identical to 
what was described with respect to Figure 5A. However, 
an electrode assembly 310 comprising two metallic probes 
320 spaced-apart perhaps 2.5 cm on a substrate 300 is 
used. Substrate 300 may be a sheet of single-sided cop- 
per clad printed circuit board measuring perhaps 5 cm x 
7.5 cm. Electrodes 320 preferably are made from brass 
and are about 0.6 cm tall, 0 . 6 cm wide, and about 1.2 cm 
in length. Spaced-apart faces of the probes define a 
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surface slanted at about 45°. Each conductive electrode 
320 is connected to one coaxial cable 12. * The finger 4 
of a subject to be. tested for glucose concentration, for 
example, is pressed against the slanted surfaces of the 
5 probes, thus completing an electrical circuit with coaxi-. 
al cables 12, and thus ports A and B of analyzer system - 
250. It is understood in Figure 5B, that port A will, 
receive back a portion of the transmitted signal as a 
reflected signal. Port B will receive that portion of 
10 the transmitted signal that propagates through the speci- 
men tissue. 

In practice, probe assembly 310 provides enhanced signal 
to noise ratio, and improved repeatability relative to 

15 other probe designs, including the probe assembly depict- 
ed in Figure 1. Reliable data have been obtained with 
probe assembly 310, typically in the frequency range of 
about 1 MHz to about 3 GHz. It will be appreciated that 
the configuration of Figure 5B is especially useful to 

20 laypersons, including suspected and actual diabetics, who 
wish to monitor their own blood chemistry, especially 
glucose concentration levels. 

Figures 6A and 6B plot predicted and actual glucose con- 

2 5 centration against time, for non- invasive ly obtained test 

data (shown by "plus signs") and for invasively obtained 
data (shown by "boxes"). Both figures depict the same 
experiment in which a human subject drank water at 14:00 
hours (2:00 P.M.) and ate food at 15:15 hours (3:15 P.M.) 

3 0 The non- invasive test data were obtained using finger 

probes 32 0 such as shown in Figure 5B, whereas invasive 
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test data were obtained from actual blood samples from 
the subject . 

Approximately 101 separate frequencies were used to ob- 
tain raw data during the experiment. Figure 6A depicts 
non-invasive predicted glucose concentration based upon ■ 
impedance and phase data taken at about 17 MHz. The 
impedance and phase data were. then converted into pre- 
dicted glucose concentration data using an algorithm. 

In Figure 6A, predicted glucose concentration shows an 
increase at about 14:20 hours, apparently corresponding 
to the Subject's intake of water. In essence, the water 
has diluted electrolyte concentration in the subject, 
which has caused predicted glucose concentration to off- 
set vertically, erroneously, by some 50 units. After 
15:15 hours, the predicted glucose level rises, which 
represents the subject's intake of food. Note, however, 
that the same 50 unit vertical offset is still present. 

Using mathematical regression analysis to examine data 
for the approximately 101 frequencies used, applicants 
realized that non-invasive phase shift data taken at 103 
MHz would provide a correction for the 50 unit error 
offset in non- invasive glucose predictions taken at 17 
MHz . 

Figure 6B shows the same experiment, now plotted with ' 
correction data taken at 103 MHz, in which "plus signs" 
depict predicted non- invasive glucose concentration data 
from the subject using 17 MHz transmission-mode impedance 
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magnitude data as corrected by the 103 MHz. phase shift 
data. Clearly the use of the higher frequency phase 
shift correction has largely compensated for the 50 unit 
offset (present in Figure 6A but not in Figure 6B) , re- 
5 suiting from water dilution of electrolytes. 

In general, Figure 6B shows close agreement between actu- 
al invasively measured glucose concentration, and non- 
invasively predicted glucose concentration. Although not 

10 fully appreciated by applicants at the time the subject 
experiment was conducted, it appears that the 103 MHz 
phase shift data provides a good measure of electrolyte 
concentration including the effects of electrolyte dilu- 
tion. At 103 MHz, small ion electrolytes including NaCl 

15 could respond to the oscillating field, whereas larger 

glucose molecules could not, and thus would not substan- 
tially influence the measurement. By contrast, the 17 
MHz data provided a measure of glucose and electrolyte 
concentration, which data when compensated for by the 103 

20 MHz electrolyte concentration cLata provided a truer mea- 
sure of predicted glucose concentration. 

Collectively, Figures 6A and 6B suggest the. wisdom of 
using data obtained at different frequencies or frequency 

25 regimes (e.g., 17 MHz and 103 MHz in this example), to 

measure different parameters (e.g., total impedance, and 
phase shift) , to provide a measure of compensation to 
more accurately arrive at the desired data (e.g., glucose 
concentration) with a greater specificity confidence 

30 level. 
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Figure 7A depicts the start lingly linear relationship ob- 
served by applicants between NaCl concentration and phase 
shift between transmitted and reflected signal at a spec- 
imen. In Figure 7A, various frequencies between 2.25 GHz 
and 2.75 GHz were used, phase difference was between two 
probes, e.g., probes 202A/B in Figure 5A. The experiment 
began with distilled water, which at shown at the bottom 
of the graph had 0 radian phase shift. Adding increments 
of 20 mM NaCl to the distilled water showed a very linear 
relationship: higher NaCl concentration increased the 
measured phase shift rather linearly. At the very top of 
the graph, data were obtained first for 300 nM NaCl, 
after which two 100 mg/dL of glucose powder was added to 
the salt water solution. As seen, in the 2.25 GHz to 
2.75 GHz frequency regime displayed, changing glucose 
concentration (indeed a rather substantial change in 
glucose concentration) did not affect phase shift mea- 
surements, whereas changing NaCl concentration produced a 
linear change in measurable phase shift. 

Figure 7B is averaged phase shift data obtained with two 
probes, using frequencies ranging from 2.0 GHz to 2.5 
GHz, in which varying concentrations of NaCl, glucose, 
and albumin were added to a baseline solution of phos- 
phate buffered saline ("PBS"). PBS was used in that it 
mimics the electrolyte environment of blood well, without 
proteins or other substances being present in the solu- 
tion. 

Consistent with the findings of Figures 6A and 6B, in- 
creasing NaCl concentration increased phase shift in a 



40 



WO 00/09996 PCT/US99/18096 
linear fashion in Figure 7B. Of special significance, 
however, is the bottommost portion of the graph, which 
corresponds to a phase shift of about 0.11 radians for a 
246 mg NaCl solution. This data line remained constant, 
5 even when 40 mg (100 mg/dl) and then 80 mg (200 mg/dl) 
glucose were added, and even when 100 mg (250 mg/dl) 
• albumin was further added. The data of Figure 7B demon- 
strates that the linear phase shift measurable for vary-, 
ing electrolyte concentration is not meaningfully influ- 
10 enced by glucose concentration and/or albumin concentra- 
tion. 

Figure 7C is a composite graph that demonstrates that a 
cross-over frequency of about 1.5 GHz renders phase shift 

15 measurements highly insensitive to varying albumin con- 
centrations in a PBS solution. In Figure 7C, the bottom- 
most-trace at about 0.017 radians represents phase shift 
caused by changing concentration of gamma globulin by 5 
g/dl, and the trace at 0.005 radians represents phase 

2 0 shift caused by changing concentration of gamma globulin 
by 2.5 g/dl. . The uppermost trace in Figure 7C represents 
phase shift due to intralipids at 1.4 g/dl concentration, 
the trace at -0.02 radians represents a different arialyte 
with glucose, not herein relevant, and the -0.015 phase 

25 shift represents intralipids at 0.7 g/dl concentration. 
Of special interest are the three tracelines centered 
about 0 radian phase shift. The trace at -0.005 radians 
represents albumin at 5 g/dl concentration, the trace at 
about -0.0025 radians represents albumin concentration of 

30 about 2.5 g/dl, and the trace at 0 phase shift is the 

baseline PBS. The various concentrations above noted are 
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substantially greater in magnitude than variations that 
would ever occur in a human being. Note that at a fre- 
quency of about 1.5 GHz , phase shift is substantially 
insensitive to albumin concentration level. Thus, by 
5 measuring different characteristics associated with a 

specimen at different frequencies or over different fre- 
quency regimes, the effects of various constituents can 
be nulled-out. In the example of Figure 7C, greater mea- 
surement specificity is attained for a desired analyte, 
10 e.g., glucose, in the presence of other substances, e.g., 
albumin. 

Figure 7D depicts phase shift data between about 300 KHz 
and 100 MHz for changing concentrations of glucose, the 
15 glucose being added to sheep blood in increments of 250 

mg%. It is seen that at about 20.1 MHz, phase shift data 
is insensitive to glucose concentration. 

Figure 8A depicts magnitude impedance data measured in a 
20 sheep blood baseline solution, for various glucose con- 
centrations, using frequencies ranging from 0.3 MHz to 
1.0 GHz. Over this extremely wide frequency regime, 
increasing concentrations of glucose increase impedance. 
The relative change of glucose concentration upon imped- 
25 ance is greater at frequencies lower than about 0.5 GHz, 
no doubt because at lower frequencies the large glucose 
molecules exert greater hinderance upon ion movement. 

In general, applicants have learned to appreciate that 
30 impedance measurement accuracy is higher at low frequen- 
cies than at higher frequencies. Thus, as will be seen, 
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impedance measurements at 2.5 GHz can provide a measure 
of glucose concentration nulling-out NaCl and other elec- 
trolyte concentrations, the equipment measurement sensi- 
tivity is substantial less than at say 100 MHz. For 
example, a measurement sensitivity of 0.1 Q is a good 
design goal. However, at 2.5 GHz, impedance magnitude • 
sensitivity will be about l/25th the sensitivity at 100 
MHz. Thus, as described herein, a recommended protocol 
will involve impedance and/or phase measurements in the 
GHz range, as well as measurements at much lower frequen- 
cies. 

Figures 8B depicts impedance change when a specimen of 
sheep's blood has glucose added, but relatively little 
change when concentrations of NaCl are added. The bot- 
tommost plot (with "boxes") is baseline sheep blood with 
a declotting agent. One addition of NaCl was then added 
(equivalent to change in concentration of 10 mM) , and 
data taken at five minute intervals for the next five 
runs. During the last (uppermost) four runs, glucose was 
added. Glucose additions clearly increase the measured 
impedance. Note that, contrary to behavior at lower fre- 
quencies, adding NaCl in the 2.94 to 3 GHz regime actu- 
ally increased impedance, probably due to an interaction 
of ions with water molecules. 

In Figure 8C, impedance data were obtaining using fre- 
quencies ranging from about 2.42 GHz to about 2.48 GHz. . 
Again, a baseline solution of sheep blood (drawn with 
"boxes") was used, into which one addition of NaCl was 
made, followed by four additions of glucose. For the 
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. NaCl additions, essentially no impedance change results 
in this frequency regime. However, the uppermost four 
runs, which represent addition of increasing concentra- 
tions of glucose, clearly increase impedance in this 
frequency regime. Thus, impedance measurements in a 
frequency regime of about 2.42 GHz to about 2.48 GHz are 
sensitive to glucose concentration, and are insensitive 
to NaCl and other small ion electrolyte concentrations. 
While sheep blood was used as the specimen, similar re- 
sults are obtainable with human blood. Further, as noted 
earlier, the human body maintains tight homeostatic con- 
trol over concentrations of most electrolytes, proteins 
and lipids within the blood. 

Figure 8D depicts impedance data for a frequency regime 
of about 2 GHz to about 2.1 GHz for a baseline of sheep 
blood (drawn with "boxes"). In the bottommost runs, the 
addition of NaCl (10 mM concentrations increments) caused 
a decrease in impedance. However, in the uppermost four 
runs, additions of glucose clearly increased impedance in 
a linear fashion. 

Figure 8E depicts impedance magnitude measurements made 
using frequencies ranging from about 2.25 GHz to 2.75 
GHz, with a specimen of distilled water into which in- 
creasing concentrations of NaCl were added. The bottom- 
most traces represent distilled water baseline data, and 
the remaining traces reflect increasing concentrations of 
NaCl, with the uppermost trace representing highest con- 
centration (200 mMNaCl). Interestingly, the effect of 
increasing NaCl concentration upon impedance varies non- 
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linearly with frequency. The right portion of Figure 8E 
demonstrates that impedance increases with increasing 
NaCl concentration (a result opposite to what is. encoun- 
tered below about 1 GHz) . By contrast, the left portion 
5 of Figure 8E shows first an increase and then a decrease 
in impedance as NaCl concentration increases (e.g., as ■ 
more Na or Cl ions are added to the test solution) . 

Figure 8F demonstrates that use of a frequency of about 
10 2.5 GHz can null -out essentially all changes in NaCl 
concentration upon impedance measurements. The data 
shown in Figure 8F were gathered using a distilled water 
specimen into which increasing concentrations of NaCl 
were added. At the approximately 2.5 GHz cross-over 
15 frequency, all curves intersected, independently of NaCl 
concentration. Note that the NaCl concentrations used in 
Figure 8F included the human physiological range of about 
135 mM to 145 mM NaCl. 

2 0 Figure 8G depicts average impedance as a function of fre- 
quency ranging from about 1 MHz to about 0 . 4 GHz Note 
that between about 0.1 and 0.2 GHz, gamma globulin ap- 
pears to saturate. 

25 In other experiments, applicants measured impedance mag- 
nitude using PBS at various temperatures to determine 
temperature sensitivity. These experiments disclosed 
that use of frequencies ' ranging from about 800 MHz to 
about 900 MHz provided impedance magnitude data that was 

30 temperature insensitive. The measurements were made 

using reflective mode, but the same result would apply to 
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transmitted mode data. When using a non- invasive in- vivo 
measurement configuration such as shown in Figure 5B, 
skin temperature at a subject's fingers can range from 
about 24°C to about 37°C. In practice, it is recommended 
that in addition to other data, that data also be taken 
in the 800 MHz to 900 MHz temperature insensitive regime,, 
to provide a measure of correction as needed for the 
other data. 

To recapitulate, the present invention recognizes that 
electrolyte ion interference, especially NaCl, with glu- 
cose measurements. can be reduced. In one application, 
the interference is effectively nulled out, using imped- 
ance magnitude measurements at a cross-over frequency. 
In another application, compensation for electrolyte ion 
effects upon glucose measurements are made. The config- 
uration of Figure 5A and likely that of Figure 5B can 
predict total glucose concentration with acceptable spec- 
ificity and error tolerance. 

As noted, it is advantageous to make high frequency and 
low frequency measurements of various parameters to pro- 
vide a good glucose concentration prediction (with good 
specificity) in a specimen. Preferably, low frequency 
regime data is taken over 21 or more frequencies, and 
high frequency regime data is taken over 81 or more fre- 
quencies. While the preferred embodiment used a network 
analyzer that provided discrete frequencies, one-at-a- 
time, the various frequencies could instead have been 
presented en masse, or as groups of frequencies, rather 
than as discrete separate frequencies. 
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High frequency, e.g., 1 GHz to perhaps 5 GHz, measure- 
ments of phase provide a good measure of electrolyte con- 
centration, in which frequency regime the phase measure- 
ments are insensitive to glucose concentration. On the 
5 other hand, use of a 2.5 GHz cross-over frequency permits 
impedance magnitude indication of glucose concentrations, 
with little contribution from electrolyte concentrations. 
But the most sensitive measures of glucose concentration 
are obtained at lower frequencies, at which impedance 
10 magnitude is a measure of glucose concentration plus 
electrolyte concentration. 

High frequency phase response was used to predict changes 
in NaCl concentration. This predicted NaCl concentration 

15 change was then used to predict the impedance magnitude 
change at low frequency due to electrolyte concentration 
change. The predicted low frequency electrolyte contri- 
bution was then subtracted from low frequency total im- 
pedance magnitude. The remainder was impedance change 

20 due to glucose concentration. In mathematical terms: 

A [NaCl] = CAL_CURVE_NaCl_PHASE_HI * APhase @ high fre- 
quency AZ NaC1 = C AL__CURVE_Na C 1 _MAG_LO * A [NaCl] 
AZ gluC08e = AZ total - AZ NaC1 
2 5 A [Glucose] = CAL_CURVEjGLUJ^AG_LO * AZ glucose 

The " CAL_CURVE " expression is derived from calibration . 
equations. When calculating concentration changes from 
phase or impedance change, it is necessary to solve an 
30 appropriate calibration equation for an unknown, e.g., 

"x" in terms of a know, e.g., "y» . NaCl calibration was 
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made using a PBS baseline solution into which NaCl was 
added in 2 mM increments, up to 12 mM above normal PBS. 
A second NaCl calibration involved diluting a PBS solu- 
tion with distilled water in 2 mM increments, to -12 mM 
from normal PBS , during which time the solution volume 
changed from about 5B8 fil to about 685 pi. However, the 
resultant calibration curve provided a linear response 
with an excellent fit, e.g., R 2 >0.999. Glucose calibra- 
tion involved three separate experiments using -10 mM 
PBS, normal PBS, and +10 mM PBS baseline solutions, into 
which glucose was added in 100 mg/dL increments to 500 ' 
mg/dL. The glucose response was quite linear with good 
correlation for the calibration curve. 

In making experimental runs, error was defined as 
100* (Predicted value - Actual value) / Actual value. On 
a run-to-run basis, NaCl concentration predictions were 
<3% and overall NaCl concentration predictions have <0.2% 
error. Overall, glucose concentration predictions had 
<13% error, and run- to -run glucose predictions had <23% 
error. These results are gratifying, although future 
embodiments will no doubt return even more accurate glu- 
cose predictions. 

The prediction method has the advantage of being fairly 
sensitive to NaCl, whose low frequency response is 
stronger than that of glucose. Although NaCl changes may 
be predicted with accuracy using high frequency phase 
data, any error in such measurement tends to be "magni- 
fied" by the leveraging effect of NaCl at low frequencies 
relative to glucose. Ideally, compensation would occur 
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at some frequency whereat the NaCl response and .glucose 
response were closer in magnitude. Applicants are also 
examining use of mathematical derivatives of the imped- 
ance and phase data obtained with the present invention. 

Modifications and variations may be made to the disclosed 
embodiments without departing from the subject and spirit 
of the invention as defined by the following claims. 
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WHAT IS CLAIMED IS: 

1. An in vivo operable method for determining 
concentration of a first chemical in the presence of a 
second substance in a specimen, the method including the 

5 following steps: 

(a) subjecting said specimen to radio frequency- 
signals having a frequency regime ranging from about 0.1 
MHz to about 5 GHz; 

(b) at a first frequency regime, using at least 

10 some of said radio frequency signals to obtain data pro- 
portional to magnitude of concentration of said second 
substance in said specimen; 

(c) at a second frequency regime, using at least 
some of said radio frequency signals to obtain data pro- 

15 portional to combined concentration in said specimen of 
said first chemical and said second substance; and 

(d) using data from said first frequency regime and 
data from said second frequency regime to obtain a mea- 
sure of concentration of said first chemical in said 

20 specimen. 

2. The method of claim 1, wherein said specimen 
includes blood. 

25 3. The method of claim 1, wherein said specimen 

includes blood, said first chemical includes glucose, and 
. said second substance includes NaCl . 

4. The method of claim 1, wherein at step (a), at 
30 least some of said frequencies are presented sequential- 
ly. 
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5. The method of claim 1, wherein at step (a) # at 
least some of said frequencies are presented 
simultaneously. 

5 6. The method of claim 1, wherein at step (b) , 

said first frequency regime ranges from about 1 GHz to : 
about 3 GHz . 

7. The method of claim 1, wherein at step (b) , 
10 said data proportional to magnitude is obtained by mea- 
suring phase shift between radio frequency signals input 
to said specimen and radio frequency signals returned 
from said specimen. 

15 8. The method of claim 1, wherein at step (c) , 

said second frequency regime ranges from about 0.11 MHz 
to about 3 GHz. 

9. The method of claim 1, wherein at step (c) , 

20 said second frequency regime ranges from about 800 MHz to 
about 900 MHz, in which regime temperature effects upon 
data are minimized. 

10. The method of claim 1, wherein at step (c) , 
25 said data proportional to combined concentration is ob- 
tained by measuring magnitude of impedance at said speci- 
men. 

11. The method of claim 1, wherein at step (d) a 
30 concentration value determined in step (b) is subtracted 

from a combined concentration determined in step (c) to 
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provide said measure of concentration of said first chem- 
ical. 

12. The method of claim 1, in which said method is 
5 carried out non-invasively on a human subject, and where- 
in step (a) includes coupling said radio frequency sig- 
nals via at least one probe that contacts a distal por- 
tion of said subject's body. 

10 13. An in vivo operable method for determining 

concentration of a first chemical in the presence of a 
second substance in a specimen, the method including the 
following steps: 

(a) subjecting said specimen to radio frequency 
15 signals at a cross-over frequency at which frequency 

concentration effects of said second substance are essen- 
tially nulled-out; and 

(b) determining from data taken at said cross-over 
frequency concentration of said first chemical. 

20 

14. The method of claim 13, wherein said specimen 
includes blood, said first chemical includes glucose, and 
said second substance includes NaCl . 

25 15. The method of claim 14, wherein said cross -over 

frequency is about 2.5 GHz. 

16. The method of claim 13, wherein at step , (b) , 
said data is impedance data. 

3 0 
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17. The method of claim 13, wherein step (a) is 
carried out non-invasively on a human subject by coupling 
said cross -over frequency via at least one probe that 
contacts a distal portion of said subject's body. 

18. An in vivo operable system for determining 
concentration of a first chemical in the presence of a 
second substance in a specimen, including: 

a transmitter outputting radio frequency signals 
having a frequency regime ranging from about 0 . 1 MHz to 
about 5 GHz; 

at least one probe [,] [coupling] coupled to said 
transmitter [, ] and adapted to contact [contacting] a 
portion of said specimen; and 

a receiver-signal processor system, coupled to said 
at least one probe, that analyzes at least some [of said] 
radio frequency signals present at said probe when said 
system is in use; 

said receiver- signal processor system providing data 
[including at least impedance and/or phase shift present 
at an interface between said specimen and said at least 
. one probe; 

wherein data provided by said receiver- signal pro- 
cessor system] that is used to determine said 
concentration of said first chemical -in said specimen. 

19. The system of claim 18, wherein said specimen 
is human blood, said first chemical is glucose, said 
second chemical includes NaCl, and wherein transmitter 
includes a network analyzer. 
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20. The system of claim 18, wherein: 
said specimen is a human subject including said 
subject's blood; 

said first chemical is glucose; 
5 said second chemical includes NaCl ; and 

said at least one probe contacts an exterior portion 
of a finger of said subject such that non-invasive data 
is provided by said system. 
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